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ABSTRACT 

o 

04 Due to its proximity, youth, and solar-like characteristics with a spectral type of K2V, e Eri is one of the most extensively 

,_H studied systems in an extrasolar planet context. Based on radial velocity, astrometry, and studies of the structure of 

^ its circumstellar debris disk, at least two planetary companion candidates to e Eri have been inferred in the literature 

^~~i (e Eri b, e Eri c). Some of these methods also hint at additional companions residing in the system. Here we present 

fvj a new adaptive optics assisted high-contrast imaging approach that takes advantage of the favourable planet spectral 
energy distribution at 4 ^ra, using narrow-band angular differential imaging to provide an improved contrast at small 

1 I and intermediate separations from the star. We use this method to search for planets at orbits intermediate between e 
r^^ Eri b (3.4 AU) and e Eri c (40 AU). The method is described in detail, and important issues related to the detectability 

^1^ of planets such as the age of e Eri and constraints from indirect measurements are discussed. The non-detection of 

companion candidates provides stringent upper limits for the masses of additional planets. Using a combination of the 
existing dynamic and imaging data, we exclude the presence of any planetary companion more massive than 3 Mjup 
anywhere in the e Eri system. Specifically, with regards to the possible residual linear radial velocity trend, we find that 
it is unlikely to correspond to a real physical companion if the system is as young as 200 Myr, whereas if it is as old as 
800 Myr, there is an allowed semi-major axis range between about 8.5 and 25 AU. 
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f^ 1. Introduction mospheric opacities cause cool objects (giant planets and 

C^ brown dwarfs) with temperatures below 800 K to emit a 

O With its distance of 3.2 pc and spectral type K2V, e Eri is i^rge fraction of their H-band flux just shortward of the CH4 

• the most nearby solar-hke smgle star. Largely due to this absorption feature at 1.6 fim, and since this feature also al- 

^ fact, It has been targeted by several attempts to mdu-ectly i^^g foj. spectral differential imaging (SDI) to enhance the 

QQ or directly detect planetary companions. So far, two exo- ^^^age contrast. However, the primary-to-secondary bright- 

Q planet candidates have been reported as companions to e ^^^^ contrast is generally smaller at longer wavelengths. 

En - e Eri b (1.55 Mjup, a = 3.4 AU), which was mferred g^,]^ ^-j^^t an improved performance can be gained for some 

from a periodic radial velocity signal by Hatzes et al. (2000) ^yp^g ^f high-contrast systems. This was demonstrated with 

and subsequently also from astrometric motion of the pri- ^n imaging survey in the L'-band of 22 young, nearby stars 

^ mary (Benedict et al. 2006), and e En c (0.1 Mjup, 40 AU), ^y Kasper et al. (2007). The brightness difference is even 

d which resuhs from an interpretation of observed irregular- smaller in M-band (see e.g. Hinz et al. 2006), but here 

ities in the dust disk surrounding the primary (QuiUen & ^^^g ^Qjgg contribution of the thermal background is also 

Thorndike, 2002). Targeted attempts have been made to significantly higher, making it observationally challenging, 

directly image both e En b (Janson et al. 2007) and e Eri jj^j.^ ^^ examine an intermediate approach, where we use 

c (Macintosh et al. 2003, Marengo et al. 2006), but did not ^ narrow-band filter towards the red end of the L'-band 

yield any detections, e Eri has also been observed as a part ^.^nge, which simultaneously allows for both good bright- 

of large-scale near-infrared adaptive optics surveys, such as ^ess contrast and observational simplicity. We test this 

Lafreniere et al. (2007), and BiUer et al. (2007). technique on the e Eri system to search for planetary com- 

TypicaUy, direct imaging surveys from the ground are panions at separations intermediate between e Eri b and e 

performed in narrow-band filters around 1.6 /im, since at- gj-j q 



Send offprint requests to: Markus Janson The paper is organized as follows: In Sect. [2] we describe 

* Based on observations collected at the European Southern the basic principles of the observing strategy chosen, and 

Observatory, Chile (ESO No. 080.C-0598). in Sect. [3] we describe the actual observations performed. 
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in terms of the techniques used and the observing condi- 
tions and corresponding instrument performance. This is 
followed in Sect. |4]by a description of the data reduction, 
and the results as given in Sect. |5.1| Our choice of filter 
is further motivated in Sect. |5.2[ where it is also explained 
over which parameter range the method should be prefer- 
able to other common methods in high-contrast imaging. In 
Sect. |5.3[ we discuss the controversial issue of the age of e 
Eri, which has fundamental implications for the detectabil- 
ity of planetary companions in the system. In Sect. |5.4| the 
brightness contrast is translated into a mass contrast. We 
examine the possible clues for an intermediate companion 
in Sect. [5T5| Finally, we summarize what these issues mean 
in terms of mass-equivalent detection limits for exoplanet 



band filter (Fnb) and its encompassing broad-band filter 
(^bb) by simply calculating: 



-2.5*logio(/NB//BB) 



(1) 



companions to e Eri in Sect. 5.6 and we conclude in Sect 

El 



2. Observing strategy 

Both the optimal choice of strategy, as well as the actual 
detection limits for a given exoplanet survey, depend cru- 
cially on the spectral energy distribution of planets in the 
mass- and age ranges of the sample. Since we lack informa- 
tion about the actual spectral distributions of young plan- 
ets from observations, we need to rely on spectral models. 
For this purpose we use the Burrows et al. (2003) mod- 
els for relative spectral distributions, and the BarafFe et al. 
(2003) models for the absolute brightness within a given 
filter. Similar strategies were applied in, e.g., Janson et al. 
(2007) and Apai et al. (2008). Here we develop these meth- 
ods somewhat further. 

The idea of choosing an optimal observing strategy from 
spectral models is based on the fact that the spectral en- 
ergy distributions of planets and low-mass brown dwarfs are 
expected to be extremely non-fiat (more precisely, deviate 
strongly from a blackbody spectrum) , with significant fiux 
windows side by side with deep absorption features in the 
near-infrared range. Optimization is achieved by choosing 
a spectral range that simultaneously contains a maximal 
amount of fiux from the planet and a minimal amount of 
fiux from the star - and, in addition, a minimal amount of 
fiux from the thermal background. Since different types of 
companions will have different spectral properties, the opti- 
mization needs to be done with consideration to what is ac- 
tually being sought after. As a general principle, a narrow- 
band filter that contains a relatively large fraction of the 
planetary fiux within the wide band that encompasses its 
spectral range will be preferable to the broad-band itself for 
the observations. Thus, for instance, the Fl narrow-band 
of SDI is vastly preferable to H-band for this purpose. This 
issue warrants emphasis, because intuition might suggest 
the opposite - for instance, in the case of an isolated faint 
object, the SNR will be limited by the amount of photons 
received from the object itself, which can be maximized by 
maximizing the spectral range of the observations. In the 
case of PSF- limited high-contrast imaging, however, the in- 
teresting signal is provided by the secondary, but the noise 
is provided by the primary. Hence increasing the spectral 
range from, e.g., Fl to H, will add little to the signal, but 
plenty to the noise, leading to a decrease of the effective 
SNR. 

Given a spectral distribution for a certain object, we can 
determine a magnitude correction Amag between a narrow- 



where /nb is the average flux density in the narrow- 
band filter, and /bb is the same in the broad-band filter. 
The fiux densities are acquired from Burrows et al. (2003). 
Under the assumption that the stellar spectral features are 
neghgible, Amag can be added to the absolute magnitude of 
the hypothetical planet under study as given by the Baraffe 
et al. (2003) models, and compared to the known absolute 
magnitude of the primary, in order to estimate the physical 
contrast corresponding to the relevant filter. The contrast 
must then be compared to the instrumental contrast that 
can be achieved corresponding to that particular filter. The 
instrumental contrast depends on a number of factors - for 
instance, due to the fact that Strehl ratio increases mono- 
tonically with wavelength, the speckle-limited instrumental 
contrast will improve with increasing wavelength. On the 
other hand, the thermal background also increases with 
wavelength, and so the background-limited instrumental 
contrast will decrease with increasing wavelength. Likewise, 
the diffraction limit of the telescope increases with increas- 
ing wavelength, which can decrease the instrumental con- 
trast at small separation. An additional difference arises 
from the fact that different techniques are applicable in 
different filters. For instance, SDI, which improves most 
PSF-related instrumental contrast is well applicable in H- 
band, but not in L'-band (at least with presently existing 
instrumentation). These effects are most efficiently taken 
into account by comparing actual measured instrumental 
contrasts, acquired under fair comparative circumstances. 
Prior to observing e Eri, we determined that NB4.05 would 
be the best filter based on the contrast curves achieved in 
L'-band (see Kasper et al. 2007), but in Sect. 5.2 we will 
use the more illustrative (and accurate) approach of de- 
riving detection limits using the actual measured contrast 
curve as presented in Sect. |5.1| 

3. Observations 

The target e Eri was observed at two neighbouring epochs, 
the first on Dec 17 and the second on Dec 28, 2007. These 
adaptive optics (AO) observations were performed in ser- 
vice mode with the NACO instrument at the VLT on 
Paranal, Chile. For the purpose of optimizing the star- 
planet contrast, the NB4.05 filter (a narrow-band filter cen- 
tered at 4.05 /im with a bandwidth of 0.02 /im) was used 
for the imaging, this choice is motivated in Sect. |5.2[ The 
pixel scale was 27 mas/pixel, which corresponds to a field 
of view of about 27" by 27". A neutral density filter which 
reduces the transmission by a factor 50 was used during the 
acquisition in order to avoid saturation in the acquisition 
image. This filter was removed during the actual observa- 
tions, which led to saturation of the detector by about a 
factor 10 at the center of the stellar PSF, and the satura- 
tion extended out to 0.2", as had been estimated before- 
hand. This does not negatively affect the result, aside from 
the fact that a few minutes need to be spent on empty 
read-outs at the end of the run in order to ensure that no 
memory effects will remain on the detector for subsequent 
observing runs. Angular differential imaging (ADI), as ap- 
plied in e.g. Janson et al. (2007), was used for this run as 
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well, since it allows for efficient PSF subtraction and yields 
easily identifiable companion signatures. This was imple- 
mented by performing the first half of the observations on 
each occasion with an instrument rotator orientation of 0° 
(with respect to North in the sky) , and the other half with 
an orientation of 33°. A jittering scheme was performed to 
allow for efficient subtraction of the strong background that 
can be expected around 4 /im. Integration times and other 
observation parameters are given in Table IT] where these 
values are also compared to the conditions of epoch 4 in the 
multi-epoch observing campaign of e Eri using NACO-SDI 
(Janson et al. 2007). 

Table 1. Observational conditions for the two runs, com- 
pared with the conditions of epoch 4 for the NACO-SDI e 
Eri campaign for reference (Janson et al. 2007). 





NB4.05, 1 


NB4.05, 


2 


SDI 


Main date 


17 Dec 2007 


28 Dec 2007 


1 Jan 2006 


Mean seeing" 


1.42" 


0.65" 




0.82" 


Mean Strehl'' 


84.3% 


87.4% 




33.8% 


Mean humidity 


6% 


32% 




38% 


Coherence time" 


2.6 ms 


3.4 ms 




3.8 ms 


Frames 


19 


19 




52 


(per angle) 










DIT 


1.0 s 


1.0 s 




1.0 s 


NDIT 


61 


61 




86 


Tot. time 


1159 s 


1159 s 




4472 s 


(per angle) 











'^ Values given by the atmospheric seeing monitor at a wave- 
length of 500 nm. 

^ Strehl ratio given by the AO system, rescaled to the observ- 
ing wavelengths. 



the whole image range. The images were also shifted so that 
the center of the stellar PSF corresponds to the center of the 
frame, where the absolute shift was determined by calculat- 
ing the flux-weighted centroid in the image and taking this 
as the saturated stellar PSF center (since e Eri is the only 
bright object in the field). Unsharp-masking was applied in 
order to remove features of low spatial frequency (stellar 
halo, background) without affecting higher-frequency fea- 
tures such as the PSF cores of actual companions. For this 
purpose, image copies convolved with a Gaussian kernel of 
FWHM 500 mas were produced and subtracted from the 
original images. Subsequently, the actual subtraction of 0° 
and 33° images was performed, followed by co-addition of 
the data from the two different nights. Since the uncertainty 
in the pixel scale is on the order of 0.1 mas/pixel (see e.g. 
Neuhauser et al. 2005), at the largest separations probed 
there can be deviations of about 50-60 mas, i.e. about half 
the FWHM. Hence, any object detected in this range could 
be expected to appear somewhat extended. 

Uncertainty levels were estimated in the same way as 
in Janson et al. (2008), by calculating the pixel- value stan- 
dard deviation in one-pixel annuli at every separation from 
the star. For translating the uncertainties into detectable 
star-planet brightness contrasts, it is necessary to estimate 
the brightness of the primary star, which can not be done 
on the basis of the science images themselves, as they are 
saturated. Instead, this was performed by estimating the 
brightness (in terms of the PSF peak value) in the acqui- 
sition frame, and rescaling to take into account the differ- 
ence due to use of a neutral density filter for the acquisition 
frame, which reduces the flux by a factor 50. The contrast 
curves of the two individual observing runs were also cal- 
culated in the same way, in order to allow the study of pos- 
sible differences in outcome due to atmospheric conditions 
or related effects. 



Although the atmospheric conditions given by the on- 
site atmospheric monitor were quite different for the two 
nights (see Table IT]) , the average quality of the wavefront 
correction, quantified by the Strehl ratio as given by the AO 
system itself, was remarkably similar. The exquisite Strehl 
ratio of well above 80 % in all cases results from the rel- 
atively long wavelengths used. For reference, observations 
performed with NACO-SDI at 1.6 /im using the same AO 
system rarely yield a Strehl ratio above 30 %. The mea- 
sured FWHM of the PSF in the non-saturated acquisition 
image, as given by the 'imexamine' routine in IRAF, was 
119 mas. 

4. Data reduction 

The data reduction was performed in a similar manner as 
discussed in Kasper et al. (2007), for the part of the sample 
presented there for which ADI was used. The basic steps of 
flat flelding, bad pixel removal and background subtraction 
were performed by the ESO 'jitter' routine in the auto- 
matic pipeline with excellent results. The remaining steps 
of ADI, combination of the images from the two different 
nights, and analysis of the data quality, were performed 
with a dedicated IDL pipeline speciflcally written for this 
purpose. For the ADI, the 0° and 33° images were shifted to 
a common position using bilinear interpolation, where the 
relative shift was determined through cross-correlation over 



5. Results and discussion 

In this section, we present the results and various consid- 
erations that go into the determination of actual detection 
limits for e Eri. First we show the actual final image, and 
corres pond ing contrast as a function of angular separation 
(Sect. 5.11. We then continue our discussion from Sect 
about the grid calculations based on theoretical spectra, 
models from the literature that were used to choose an op- 
timal filter for the observations, and how (in combination 
with the measured contrast) they yield detection limits for 
the whole range of ages and masses allowed by the models 
(Sect 



5.2 1. Since the age of the system is fundamental in 



constraining these detection limits, an attempt to constrain 
the age as fa r as possible based on the existing literature 
follows (Sect. [5!3| . The brightnes s con trast is subsequently 
translated into mass limits (Sect. 5.4 1. Further constraints 
are then added from independent dynamical data (Sect. 
5.5). Finally, all of these issues are combined to produce 
detection limits and constrain the allowed paramete r ran ge 
of any planets that may reside in the system (Sect. 5.6). 



5.1. Difference image and contrast 

The final output image is shown in Fig. [T] The remain- 
ing residuals in the image appear to be primarily due to 
rotation of PSF substructures between 0° and 33°, and be- 
tween the first and second observing run. Ideally, all the 
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individual observations should be carried out at identical 
or symmetrical telescope orientations, for instance through 
symmetrical distribution of on-target time around the time 
of minimum zenith angle during a night. SDI data tends 
to efficiently mitigate effects such as spider rotation due 
precisely to the fact that the spectral image pairs are, by 
definition, taken at identical orientations. In any case, the 
achieved image contrast is excellent, as shown in Fig. [2] For 
instance, at an angular separation of 2.6" we reach a con- 
trast of A7714.05 = 13.2 mag, which corresponds to a star- 
to-planet flux ratio of 1.6*10^. The contrast corresponds to 
the detectability of a companion by 3 a. While this is not a 
very statistically stringent detection constraint for a single 
putative signature in the innermost part of the image where 
systematic noise dominates (see Fig. |3 1 , it is stringent for 
the vast majority of the field where the residual noise is 
uncorrelated. The ADI technique in our application also 
yields two independent signatures, which increases the rel- 
evance of a putative detection. Thus we follow the previous 
convention for this technique and use 3 cr as the criterion. 
The difference in performance between when ADI is used 
and when it is not (for a single object signature) is shown 
in Fig. |4] 




4 6 8 10 12 

Angular separation (arcsec) 

Fig. 2. Image contrast achieved in the NB4.05 filter with 
the ADI technique. Solid line: The full combined data set. 
Dashed and dotted lines: The two individual observing runs 
before co-addition. The decrease in contrast beyond about 
10" is due to the fact that the individual jittered frames 
overlap completely out to 10", but non-completely out- 
wards of this separation. Hence, the available effective ex- 
posure time decreases with separation in this range. 



Also shown in Fig. [2] are the contrasts of the individual 
observing runs. The results from the individual runs are 
very similar, which is consistent with the similar Strehl ra- 
tios in the two cases. The contrast for the full dataset is sig- 
nificantly better than for the individual ones, implying that 
the residual noise is mostly incoherent, and that we benefit 
from using all the available data, which is not inherently ob- 
vious when the atmospheric conditions are so different. At 
4.05 /im, the isoplanatic angle is typically larger than 30", 
hence anisoplanatic effects are not expected to significantly 
affect the quality anywhere in the image range. 




Fig. 3. Six artificial companions at 3(t significance, two 
each at separations of 1", 2", and 3". The center of the 
image has been artificially set to zero to enhance visibil- 
ity. 3cr is a useful constraint for most of the image range. 
However, in some cases it is problematic, such as for the 
rightmost dark signature, which by chance coincides with 
a local maximum in the residual stellar PSF and becomes 
difficult to distinguish. North is up and East is to the left 
in the image. 



o 
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Fig. 4. The difference in performance between when ADI is 
used and when it is not. A positive number means that the 
contrast has increased when using ADI, for a single object 
signature. There is an apparent decrease in performance at 
large separations - this is because there was no systematic 
noise in that range to be subtracted out, whereas random 
noise increased in the (negative) co-addition performed dur- 
ing ADI. However, at the same time, the number of object 
signatures is doubled, so there is in reality no decrease of 
performance in this range. The most substantial increase in 
performance is, as expected, at small angular separations, 
where systematic noise sources dominate. 



5.2. Comparison with other techniques 

Detection limits for different methods in terms of planetary 
mass are derived in the following way: For each sampled 
angular separation and each filter, the measured contrast 
is compared with a physical contrast from the Baraffe et 
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Fig. 1. Output difference image of e Eri. The orthogonal line-shaped artifacts toward the edges of the image are due 
to the jittering, which shifts the field of view by a few arcseconds for each individual integration. In an unsaturated, 
undifFerenced image, the PSF peak of e Eri would here correspond to 3.5 * 10^ counts at a Strehl ratio of 85 %. North is 
up and East is to the left in the image. 



al. (2003) models, as described above. The comparison is 
made through linear interpolation between the grid points 
of the models to keep the result smooth - for instance, 
if a 5 Mjup planet requires a contrast of 12 mag, and a 
6 Aijup planet requires 11 mag, a measured instrumental 
contrast of 11.5 mag yields a detectable mass of 5.5 Mjup. 
In the case of broad-band filters, this gives immediately the 
detectable mass at that separation. For narrow-band filters, 
we additionally need to apply the Amag correction using the 
Burrows et al. (2003) model. This has to be done iteratively, 
because a certain input mass gives a certain Amag, which in 
turn gives a new mass, etc. After 4-5 iterations, convergence 
is reached and the mass limit thus acquired. 

For the case of SDI, we applied an additional correc- 
tion. When qualitatively describing the process of SDI, it is 
often assumed that a planet with CH4 absorption contains 
no flux within the ^3 filter. In this way, when performing 
the Fl — F3 subtraction, it is easily understood that the 
stellar flux is canceled out, whereas the planetary flux is 
unaffected by the operation. This is, of course, a simplifi- 
cation, and not even a good approximation in many cases, 
such as for young and massive planets that do have sig- 
nificant flux in F3. One might think it more appropriate 
to regard the remnant flux as the difference between the 
individual Fl and F3 fluxes, but this is also an excessive 
simplification. What really happens during SDI is that the 
full Fl and F3 images are rescaled to a common A/D-scale 
in order to optimize the subtraction of the primary PSF. 
As a part of this process, the F3 position of the prospec- 
tive planet is shifted with respect to the F3 position by 



a factor A3/A1 = 1.625/1.575 « 1.03. Hence, for a given 
angular separation p between primary and secondary, the 
relative shift will be 1.03p - p = 0.03p, so if 0.03p < X/D, 
then most of the -F3 fiux from the hypothetical companion 
will be subtracted from its Fl flux. If, on the other hand 
0.03/9 ^ X/D, then F3 will yield an independent signature, 
such that is adds to the achievable contrast. Around A = 1.6 
/im, the turnover point is on the order of pto ~ 1-3" for the 
VLT. As implied, the physical reason for this is that beyond 
Pto, there will be two independent signatures (one of posi- 
tive counts, one of negative) of the companion in the final 
SDI image, separated by more than the FWHM of the PSF 
core. Hence for this special case, it is easier to detect those 
companions that do not exhibit significant CH4 absorption 
than those that do. In our simulations, we take this effect 
into account, considering both its dependence on angular 
separation and spectral energy distribution, with a correc- 
tion factor applied to the contrast curve. The correction 
is calculated by generating two Gaussian (ID) functions 
with a FWHM set by X/D, a relative shift given by 0.03p, 
and a relative amplitude fiuxes in the Fl and F3 filters, 
respectively. The Gaussians are then subtracted, and the 
correction factor is calculated as the sum in quadrature of 
the negative and positive residuals over the absolute sum of 
the Gaussian corresponding to Fl. A^ag between Fl and 
H is then calculated as described above. 

Each of the calculations are made corresponding to a 
range of ages. The Burrows et al. (2003) models exist down 
to 1 Mjup for the ages of 100 Myr, 320 Myr and 1 Gyr. 
For these ages, the models are complete up to 10 Mjup. 
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Hence, we set the boundaries of our simulations to 0.1-1 
Gyr and 1-10 Mjup, since the BarafFe et al. (2003) models 
are also complete in this range. By complete, we mean that 
grid points exist at or around these boundaries such that 
no extrapolations need to be made. 



Performing all the calculations described above, for all 
relevant masses, separations, ages, methods, and wave- 
length ranges, we can calculate detection limits and com- 
pare different methods. While we have looked at several dif- 
ferent filters available with NACO (e.g. broad-bands such as 
H and Ks, and narrow-bands such as IB2.09 and NB2.17), 
here we focus on the three most interesting cases: NACO- 
SDI, L'-band and NB 4.05. In Fig. [5] we show a compar- 
ison between ADI in NB4.05 with NACO, and SDI+ADI 
around 1.6 //m with NACO-SDI, for the ages 100 Myr, 320 
Myr, and 1 Gyr. Since the instrumental contrast curve for 
NB4.05 is jagged due to diffraction ring residuals, the mass 
limits that are based on this contrast become jagged as well, 
which makes the diagram more difficult to read out. For this 
reason, we have neighbour-averaged each separation point 
with the two nearest neighbours in each direction. This is a 
purely cosmetic operation which does not affect the factual 
content of the figures. Also shown is ADI corresponding to 
L'-band with NACO. In that case, the contrast is not based 
on a direct measurement, but on the assumption that the 
NACO ADI performance is equal in NB4.05 and L'-band 
(i.e., any difference is entirely due to the spectral energy 
distribution of a given object, which is generally favourable 
in NB4.05). This is a fair assumption, since the method and 
instrument is the same, though minor differences will exist 
since NB4.05 is not at the center of the L'-band, and thus 
the average Strehl ratio, background, and diffraction limit 
will be slightly different. The most important consequence 
of this result is that if the theoretical models are to be 
trusted, then ADI in NB4.05 significantly outperforms all 
other tested methods for all relevant masses and ages. This 
statement holds as long as the target has a similar near- 
infrared brightness as e Eri. If the target is significantly 
fainter, the background will start to dominate, and eventu- 
ally the contrast gained with respect to, e.g., NACO-SDI 
(which is not significantly affected by the thermal back- 
ground at this level), will be void and the latter will be 
preferable. Another point that bears mentioning is that the 
results are specific to the AO system of NACO, and the 
Strehl ratio is one of the most important parameters gov- 
erning the relative results at 1.60 and 4.05 /im. With a 
next-generation AO system, the Strehl ratio improvement 
at 1.60 /xm can be expected to be large, whereas the Strehl 
ratio improvement at 4.05 /zm, where the ratio is already 
very high, will be minor. For this reason, it is possible that 
SDI+ADI at 1.60 ^m might still outperform ADI at 4.05 
/xm for the next-generation planet-finding instruments. 



As a final point concerning the relative performance be- 
tween NACO at 4.05 /im and NACO-SDI, we note that the 
field of view is considerable larger in the former case, and 
thus additionally better suited for the study of more distant 
companions. In the comparative plots, we have cut off the 
outer range in order to highlight the inner range, but note 
that NACO at 4.05 /im extends out to 10" and beyond (see 
Fig.[2]). 



5.3. Age of e Eri 

As we could see in Sect. |5.2| the detectable mass is strongly 
dependent on the age of the system. Determining the age of 
young stars is however very difficult, and this is in particular 
true for the case of e Eri. In fact, age estimations in the 
literature have been as disparate as the extent of our grid, 
about 0.1-1 Gyr. Here we will briefly examine this issue, 
to see if we can say something about the relative merits of 
different age estimations, and how far we can constrain the 
age range. 

The methods that have been used to infer an age of 
e Eri are large in amount: Isochronal dating, kinematics, 
Ca II activity. X-ray activity, Li abundance, and rotation. 
Unfortunately, these methods give widely discrepant re- 
sults, and many, if not all, have indications of being not 
well applicable for the particular case of e Eri. In our dis- 
cussion, we will focus on three papers that concern dating 
of e Eri, Song et al. (2000), Fuhrmann (2004), and Barnes 
(2007), that cover all these methods. While several other 
papers exist that address the issue (e.g. Lachaume et al. 
1999, Di Folco et al. 2004, Saffe et al. 2005), they cover 
the same methods and reach similar conclusions for each 
method, and the three former papers appear to be the most 
appropriate for age determination in this context. 

Fuhrmann (2004) uses kinematics to associate e Eri with 
the Ursa Major association, giving an age estimate of 200 
Myr. Barnes (2007) uses gyrochronology (stellar rotation) 
to get an age estimate of 440 Myr. Song et al. (2000) 
ambitiously use five different methods (color-magnitude- 
diagram (CMD) with isochrones, Li abundance. X-ray 
emission, Ca II chromospheric activity, and young/old disk 
kinematics) and quote a final age range of 530 to 930 Myr. 
At a first glance, it may appear that the Song et al. (2000) 
results are far more secure, with five different methods con- 
sistently yielding an older age than the two other papers. 
However, we need to carefully check the applicability of 
each of the different methods. First of all, as is clearly noted 
also by Song et al. (2000), isochronal dating is not applica- 
ble to a target like e Eri, because there is so little difference 
between different ages in this spectral type range beyond a 
few tens of Myrs. For instance, it is clear from the Song et 
al. (2000) CMD analysis that the age from this method is 
not only consistent with ZAMS, but it is also entirely con- 
sistent with the Pleiades empirical isochrone. Hence, the 
isochronal method does not constrain the age and is thus 
not useful for the analysis. The young/old disk kinematics 
analysis in the same paper identifies e Eri as an old disk 
star, but it is very close to the border that separates the 
old and young disk, and the uncertainty of the identifica- 
tion procedure is not shown. Since the kinematic analysis 
of Fuhrmann (2004) identifies the spatial motion of e Eri as 
being consistent with the Ursa Major association, and the 
analysis of that paper is more ambitious in this particular 
sense, we conclude that the kinematics of e Eri does not ex- 
clude a younger age than the old disk population in general, 
and disregard this part of the Song et al. (2000) analysis. 
What remains are X-ray emission, Ca II activity and Li 
abundance, all of which are consistent with an age similar 
to the Hyades of 600-800 Myr. However, problems remain. 
Israelian et al. (2004) show that stars in the e Eri spec- 
tral type range with detected exoplanets are significantly 
Li depleted compared to a reference sample. This result is 
confirmed by Gonzalez (2008), who finds the same trend 
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also for rotation and activity. This implies that the val- 
ues of Song et al. (2000), and even of Barnes (2007) should 
rather be seen as upper limits than realistic age estimations. 
Furthermore, there is a comment to be made with regards 
to the discrepancy between the age estimations from rota- 
tion (440 Myr) and activity (600-800 Myr). The physical 
link between activity as shown by Ca II and X-ray data 
is not well understood (Barnes 2007), but the activity is 
linked with the magnetic dynamo, which is in part gov- 
erned by the rotation rate of the star. In this context, there 
is at least one obvious connection between age and activity: 
Stars are known to lose angular momentum with age, hence 
once they have reached the main sequence, the rotation rate 
decreases, the dynamo strength follows, and the activity de- 
creases as a consequence. Note that in such a framework, 
activity is a less direct age indicator than rotation rate. 
This makes the former less reliable than the latter. This is 
more than a hypothetical line of reasoning, since it is known 
that activity cycles exist in solar-type stars on timescales 
over which the rotation rate can be approximated as being 
constant. Such cycles even include extended periods of in- 
activity - for instance, during the period ~1645-1715, the 
Sun showed almost no sunspots (the so-called 'Maunder 
minimum', after e.g. Maunder 1922). Meanwhile, if the ro- 
tation rate develops in a non-standard way, then within the 
same framework, both indicators would be unreliable. 

Kinematic group membership as applied in Fuhrmann 
(2004) is a method that differs qualitatively from the other 
methods mentioned. If a given star is correctly identified 
as a member of a group, the age estimate is more precise 
than for any of the other methods, but on the other hand, 
if the star is misidentified and shares a common proper 
motion only by chance, it is irrelevant for dating purposes. 
The identification of e Eri as an Ursa Major association 
member is not very convincing, as it is one of the most 
deviant stars in terms of spatial motion compared to the 
mean motion of those considered as group members. Thus, 
while we can't rule out e Eri as an Ursa major association 
member, we also don't consider it a necessary match, and 
so we only use 200 Myr as a lower limit to the age range. 

In summary, following the reasoning above, we consider 
the most probable age of e Eri to be around 440 Myr (i.e. 
the age estimated from rotation rate in Barnes 2007), but 
we will use an age range of 200-800 Myr which corresponds 
to all the methods that we consider to be of some relevance 
for the age estimation. 



5.4. Mass limits 



Using the simulations described in Sect. 5.2 and the age 
estimate from Sect. |5.3| we get mass-equivalent detection 
limits as a function of physical separation shown in Fig. [6] 
The radial velocity trend for which a corresponding curve 
is also shown in that figure is discussed in Sect . |5.5| and the 
consequences for detection limits that arise from this will 
be discussed in Sect. 15.61 



5.5. Input from dynamical methods 

Direct imaging data can be combined with data from dy- 
namical methods to constrain the presence of planetary 
mass companions. Radial velocity and astrometry data have 
already been used to infer the presence of e Eri b (Hatzes 



et al. 2000 and Benedict et al. 2006). This input was used 
for the images in Janson et al. (2007) and is not a con- 
cern of this paper (the predicted projected separation is 
too small at the epoch of the NB4.05 images to improve 
on the SDI contrasts, see Sect. 5.6 1 ~ here, we are instead 
primarily interested in what could be said about additional 
companions. In addition to these conventional methods, e 
Eri offers an additional option to study dynamics in the 
system through its spatially resolved circumstellar debris 
disk, which can be expected to be shaped by interactions 
with any massive planets residing in the system. 



5.5.1. Debris disk structure 

Using numerical modeling, Quillen & Thorndike (2002) find 
that the observed properties of the dust disk can be ex- 
plained through the presence of a planetary body with ap- 
proximate properties Mpi = 0.1 Mjup, Opi — 41.6 AU, and 
gpi = 0.3. In a more detailed set of simulations, Deller & 
Maddison (2005) confirm that this is a feasible explana- 
tion, but claim that an additional, more massive planet is 
required further in to explain the inner cavity of the disk. 
This additional planet is not simulated in sufficient detail 
to give specific predictions, but the authors give a prelim- 
inary parameter range of 10-18 AU and a mass around 1 
Mjup. Hence, the result could be indicative of an additional 
planet candidate intermediate between the two other can- 
didates (e Eri b and e Eri c). To simplify further discussion 
but simultaneously emphasize the speculative nature of this 
candidate, we will refer to it as e Eri x. 



5.5.2. Radial velocity 

In the best RV fit of e Eri b, there is a residual linear trend 
over the ^^27 year baseline, possibly indicative of a long- 
period physical companion. The trend is too large to be 
explained by secular acceleration (Benedict et al. 2006). It 
is also far too large to be caused by e Eri c - at a semi-major 
axis of ^ 40 AU, the orbital period would be ~ 280 years, 
and the resulting amplitude of the curve would correspond 
to a planetary mass of ~ 8 Mjup, assuming that the orbital 
plane of the planet is the same as the plane of the disk, 
the rotational plane of the star, and the predicted orbital 
plane of e Eri b, i.e. has i « 30 ° (Saar & Osten 1997, 
Greaves et al. 1998, Benedict et al. 2006). This is almost two 
orders of magnitude in mass above the Quillen & Thorndike 
(2002) prediction, and should have been readily detected by 
previous imaging surveys (see Marengo et al. 2006). For e 
Eri X, the inner possible range around 10 AU corresponds 
to a very short orbital period (35 years) with respect to the 
baseline of the linear trend (27 years) , but the outer range 
of 15-18 AU yields reasonable periods of 65-85 years, and 
masses around 1-2 Mjup, i.e. entirely consistent with the 
tentative values in Deller & Maddison (2005). 

5.5.3. Astrometry 

In Wielen et al. (1999), e Eri was classified as a so-called 
A/i binary candidate at the limit of detectability (roughly 
3(t). A/i binarity means that the short-term proper mo- 
tion as measured by Hipparcos is different from the long- 
term proper motion given by the catalogues FK5 (the fifth 
Fundamental Catalogue) and GC (the General Catalogue), 
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or a combination of the FK5 and GC positions with the 
Hipparcos position. The presence of a physical companion 
is a hkely cause of such a trend, because an apparent short- 
term proper motion will include the astrometric effect of a 
sufficiently massive unseen companion, whereas over a long 
baseline of up to about 200 years such as in the FK5 case, 
a companion with an orbital period of a few decades or 
less will average out, and the measured long-term proper 
motion will correspond more strongly to the motion of the 
system as a whole. Since there is now a new data reduction 
made for Hipparcos data which significantly decreases as- 
trometric residuals for bright stars like e Eri (see e.g. van 
Leeuwen 2007), we have repeated the Wielen et al. (1999) 
analysis for this target with the new improved data. 

We summarize the A/x quantities along with the corre- 
sponding significance in Table [2] The reader is referred to 
Wielen et al. (1999) for specific definitions and discussion of 
the quantities involved - we just briefly summarize the con- 
cepts relevant to our discussion in the caption of Table |2] 
For an intermediate orbit (P ~ 30 yr) as we are interested 
in here, Fqf should be small, and the other values should 
be large in order to indicate (with a corresponding level 
of confidence) the presence of a companion. Note that /io 
and /io(GC) did not change when using the new reduction 
of the Hipparcos catalogue, since the positional difference 
between the old and new reduction is so small that it does 
not make a difference for the derived proper motion. 

Table 2. A/u binary parameters with the old and new 
Hipparcos data reduction. 



Quantity 



Old red. 



New red. 



A/iFH,Q 


[mas/yr] 


1.70 ±1.08 


0.43 ±0.50 


A/iFH,5 


[mas/yr] 


3.31 ±1.05 


1.80 ±0.56 


Ffh 




3.26 


3.40 


A/iOH.a 


[mas/yr] 


1.76 ±1.02 


0.49 ±0.37 


A/iOH,i 


mas/yr] 


2.49 ±0.96 


0.98 ±0.35 


FoH 




2.78 


3.28 


A/iOF.a 


mas/yr] 


0.06 ±0.54 


0.06 ±0.54 


A/iOF,i 


mas/yr] 


-0.82 ±0.60 


-0.82 ±0.60 


FoY 




1.39 


1.39 


A/io{GC) 


H,a [mas/yr] 


1.40 ±1.03 


0.13 ±0.39 


A/io{GC) 


H,5 [mas/yr] 


3.82 ±0.96 


2.31 ±0.38 


-Fo(GC)H 




3.96 


6.22 



Note: A^ is the difference between diflerent proper mo- 
tion values identified by the indices, where 'F' denotes FK5 
proper motion, 'H' denotes Hipparcos proper motion, and '0' 
the change in position over baseline for FK5 and Hipparcos. 
0(GC) is the same as the latter, but for GC instead of FK5. 
fpH, -FoH, FoY and -Fo(GC)h are the corresponding signifi- 
cance parameters, where _Ffh = 3.44 corresponds to a 99.7 
% significance. 



It can be seen from Table |2] that all the relevant signifi- 
cances, which were already high enough for Wielen et al. to 
consider e Eri a limit-case A/i binary, have increased with 



the new reduction. In addition, Ft 



0(GC)H 



6.22 corresponds 



to a false alarm probability less than 4 * 10~^, indicating 
a strong detection of a companion. However, there are rea- 
sons to treat these numbers with great caution. In particu- 
lar, note that in addition to decreased errors, the A/x values 



have also decreased drastically - this is why the FK5-based 
significance only increases marginally. It is also surprising 
that the A/ion and the A/io(GC)H are so different, since they 
should correspond to roughly the same observable, just for 
different long-term catalogs. For these reasons, we conclude 
that while there are indications from Hipparcos for a phys- 
ical companion that could be compatible with the e Eri x 
hypothesis, we emphasize that it is far from conclusive. 

5.6. Detection limits 

In Fig. |6] we plot the expected mass of the hypothetical e 
Eri X, under the assumption that it is responsible for the 
linear RV trend, and is on a circular orbit co-inclined with 
the rest of the system (including e Eri b) . Under these as- 
sumptions, we can exclude the e Eri x hypothesis from 14 
AU and outwards if the age is 200 Myr, and from 22.5 AU 
and outwards if the age is 800 Myr. Including also the con- 
servative constraint that the period has to be substantially 
larger than twice the baseline, we find that if the system 
is 200 Myr old, e Eri x cannot exist, whereas if it is older, 
there is an allowed parameter range between 14 and 22.5 
AU, and between 1.0 and 2.7 -Mjup. 

Unfortunately, exoplanet observations have shown that 
a perfectly circular orbit is not an ideal approximation - 
many of the detected exoplanets have significantly eccen- 
tric orbits. Such orbits complicate the analysis due to their 
time-dependency in physical separation and orbital veloc- 
ity. A planet on an elliptical orbit near periastron will reside 
further in towards the star for a given semi-major axis, and 
its mass corresponding to a given radial velocity amplitude 
will be smaller compared to a planet at the same physi- 
cal separation and orbital phase, but with a circular orbit. 
This is due to the fact that the planet near periastron will 
move faster, such that a smaller mass is required to yield 
the same amplitude of the stellar reflex motion. However, a 
planet observed at a non-special (i.e., random) epoch will be 
unlikely to reside near periastron, because it spends a very 
small fraction of its orbit in that phase. In addition, the RV 
trend will change on a much faster timescale for a planet 
near periastron, such that to explain a linear trend with a 
baseline of 27 years, the full orbital period has to be several 
times larger, requiring a large orbit and thus also a large 
planetary mass. For a planet near apastron, the opposite is 
true. A planet on an elliptical orbit is likely to reside near 
apastron at a random epoch, and for a very eccentric orbit, 
the orbital period only needs to be a rather small amount 
(< 2 times) larger than the observed baseline. On the other 
hand, the lower orbital velocity compared to a planet on a 
circular orbit with the same physical separation leads to a 
larger mass to yield the same velocity amplitude. 

To estimate the impact of orbital eccentricity on our 
detection limits, we perform a simple simulation in the 
following way: We wish to simulate detection probability 
Xdot as a function of semi-major axis a, and for this rea- 
son we step through 0.1 to 30 AU with steps of 0.1 AU. 
A representative eccentricity distribution is chosen by us- 
ing the observed eccentricity distribution of RV-detected 
exoplanets (from th e Extrasolar Planets Encyclopaedia, 
http://exoplanet.eul outside of 0.1 AU where orbital cir- 
cularization ceases to be important. This gives 179 realistic 
eccentricities. For each semi-major axis and eccentricity, we 
simulate 200 uniformly random orbits and instances, rep- 
resented by the argument of periastron uj and the fraction 
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of the period /p — t/P, where < t < P. Following the 
standard equations for Keplerian motion, this yields the 
eccentric anomaly E (since the dependence £'(/p) is non- 
linear, this is done numerically), which in turn yields the 
projected separation p. Xdot is determined as 



Xr 



if Af,, 



{ai)/wi > Mc. 
otherwise 



ipd 



(2) 



where i is an index representing each simulated orbital 
instance, Mcdpt) is the measured detection limit at the 
calculated projected separation and Mi-v(ai) is the mass 
predicted by the RV trend for the a at the given i if the 
orbit was circular. Wi = Vi/vo,i is the fraction of real orbital 
velocity at the epoch corresponding to i over the velocity 
of a circular orbit with semi-major axis a^ - i.e., it properly 
takes into account the eccentricity of the orbit. In addition 
to this, we must ensure that simulation i is consistent with 
the baseline of observations. This is difficult to do strin- 
gently, as one would need to have specific information about 
the curvature of the (quasi-) linear trend, which would be to 
over-interpret the data at hand. In order to remain flexible 
with regards to curvature, we set the requirements that: 1) 
P{ai) > 2 * 27 * Wi * q yr, and 2) P{ai) > 27 * q yr, where 
q is an arbitrary parameter governing which fraction of a 
period the baseline can correspond to. When q = 1, the 
first condition claims that for a circular orbit, the baseline 
covers at least half the period (which is clearly as conserva- 
tive as it gets), and that if the eccentricity is non-zero and 
the companion is close to periastron (wi > 1), then the pe- 
riod must be larger because the observed trend represents 
a smaller fraction of the observed orbit, and the opposite 
holds true for apastron. The second condition at q = 1 ac- 
counts for the fact that no matter how close to apastron a 
companion is, it clearly could never have a period shorter 
than the baseline itself (again, as conservative as it gets). 
By setting q to some value g > 1, we can adjust our limits 
to be less conservative and, to some extent, more realistic. 
Here we will consider two different cases: g = 1 and q = 2. 

Finally, we can calculate Xdot as a function of tti as the 
mean of all Xi corresponding to a given at. Only the Xi 
that pass the period requirements are included in the mean. 
Separations a^ that are small enough that no Xi pass the 
criterion, such that no mean can be evaluated, are not in- 
cluded. We plot Xdet for the ages of 200 Myr and 800 Myr, 
and g- values of 1 and 2 in Figs. [7] and [8] We see that there is 
now an allowed parameter range between 8.5 and 16.3 AU 
(for Xdct < 0.95) for e Eri x even at 200 Myr for q = 1, this is 
due to small orbits near apastron. We note that this param- 
eter space is almost entirely due to the very conservative 
(or, from the perspective of allowed near-apastron orbits, 
very generous) period limits. If a more intuitive constraint 
such as 5 = 2 is set, then Xdet = 1 everywhere. For 800 
Myr, the corresponding ranges are 8.5 to 27.3 AU for g = 1 
and 8.5 to 24.7 AU for q = 2. 

We note that regardless of whether there exists an e 
Eri X or not, the RV trend serves as a useful indicator of 
what sort of companions could be detected in the RV data. 
For instance, returning to Fig. [6] we see that at 8.5 AU 
where the 800 Myr contrast curve starts to be sensitive to 
3 Mjup planets, the circular RV trend would correspond to 
less than 0.5 Afjup, so if there was a 3 Afj^p in the sys- 
tem at that separation, there would be an RV trend with 
a velocity amplitude 6 times larger than the existing linear 



trend, which would make the amplitude even larger than 
the one attributed to e Eri b, and extend over a much longer 
baseline, making it very easy to discern. For smaller semi- 
major axes, the amplitude would be even larger. Thus since 
no such trend exists, we can combine the RV and imaging 
data to conclude that there are no companions to e Eri 
more massive than 3 Mjup. This is valid over all possible 
semi-major axes, since we are sensitive to 3 Afjup out to 14", 
outward of which Marengo et al. (2006) provide a sensitiv- 
ity to 1 Afjup companions. Fig. [6] corresponds to a circular 
orbit, but the conclusion can be extended to elliptical or- 
bits, since the only way a planet can be more massive and 
still correspond to the same RV amplitude at a given semi- 
major axis is if it is near apastron, in which case it would 
be easier to detect in the imaging. 

Of course, there is a hypothetical possibility that a com- 
panion could exist that does not share an orbital plane with 
the rest of the system - in that case, projection effects could 
still mask a massive companion in rare cases. However, for 
the case of disk-formed planets, which we would certainly 
expect to share an orbital plane with the circumstellar disk, 
we can conclude that whatever the actual number and or- 
bital parameters of planetary companions to e Eri are, they 
all must be less massive than 3 Afjup. 



6. Conclusions 

We have performed a deep, high-contrast imaging observa- 
tion of e Eri with a new method that allows for detection of 
lower-mass planets at small separations from bright stars 
than previously possible. The method is based on ADI in 
the NB 4.05 /im filter, and compares favourably to other 
methods in the literature. Its advantage is however lim- 
ited to very bright stars, since the background limit will 
dominate the contrast limit for faint stars, such that other 
methods (e.g. SDI+ADI or L'-band ADI) are preferable for 
such targets. As the final output images show no signs of 
any convincing companion candidate, we place constraints 
on any planets that may reside in the system. In partic- 
ular, we test the hypothesis of the existence of an inter- 
mediate orbit planet possibly indicated by radial velocity, 
Hipparcos A/i astrometry and numerical attempts to re- 
produce the circumstellar debris disk structure of e Eri. In 
order to emphasize the hypothetical nature of such a com- 
panion, it is referred to as e Eri x. We find that, due to 
the uncertain age of the e Eri system, which we quote as 
200-800 Myr based on a detailed investigation of the exist- 
ing literature, the presence of e Eri x can not be excluded 
based on our images, but its allowed parameter range can 
be constrained. In general, a combination of the existing 
RV data and high-contrast images allows us to conclude 
that whatever the amount and distribution of planets is in 
the e Eri system, none of them can be more massive than 
3 Mjup. 
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Fig. 5. Comparison of different high-contrast methods in 
terms of detectable mass of companions (in the simulated 
range of 1-10 Afjup) to stars of ages 100 Myr (upper panel), 
320 Myr (middle panel), and 1 Gyr (lower panel), and of 
similar brightness to e Eri. Solid line: 4.05 /xm with ADI. 
Dashed line: NACO-SDI with ADI. Dotted line: L'-band 
(assumed, see text) with ADI. Note that the NACO-SDI 
curve had about twice the integration time as the others. 
Under the assumption that SNR ~ i^/^ (see Janson et al. 
2007), this means that under the same integration time, 
the relative performance of NACO-SDI would get worse by 
0.36 mag, generally corresponding to < 0.5 Aijup in these 
cases. Also note that L'-band is probably preferable to NB 
4.05 in the background-limited regime, but the background 
limit is not yet reached at the outer cut-off of the figure 
(compare with Fig.l2|. 
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Fig. 6. Detection limits for ages of (from top to bottom) 
800 Myr, 600 Myr, 400 Myr, and 200 Myr. Also plotted is 
the mass corresponding to the linear RV trend (solid line 
increasing with separation) given the e Eri x hypothesis as 
described in Sect. |5.5| and a vertical dotted hne correspond- 
ing to the measured baseline of the trend. The dashed line 
is the minimum projected separation ai i — 30° . The latter 
curve serves to show that the detection limits are flat over 
the projection range for all ages, hence the projection is ir- 
relevant and e Eri x would simply be detectable at a given 
age from the separation where the solid RV line crosses the 
relevant detection limit, and outwards. 
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Fig. 8. Same as Fig. [t] but at an age of 800 Myr. 
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Fig. 7. Detection probability as a function of semi-major 
axis, at an age of 200 Myr. Thick lines: Xdet for 5 = 1 (solid 
line) and q = 2 (dashed line). Thin lines: Fraction of oc- 
currences for a given semi-major axis that yield acceptable 
orbits, in the sense that the hypothetical planet could be 
consistent with the basic features of the RV trend (such as 
minimum baseline) , for q — 1 (solid line) and q = 2 (dashed 
line). The meaning of the latter curves is that in the ranges 
where they give low probability values, it is unlikely that a 
planet could reside at those semi-major axes and still have 
an orbit that is consistent with the measured RV trend, 
according to the conditions specified by the q values. 



